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Abstract. Shapes of relative emission sources can be accessed by expanding 
shapes of correlations at low relative velocities in pair center of mass in Carte- 
sian harmonics. Coefficients of expansion for correlations are related to the 
respective coefficients of expansion for the sources through one dimensional 
^ , integral transforms involving properties of pair relative wavefunctions. The 

■ methodology is illustrated with analyses of NA49 and PHENIX correlation 

data. 
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1. Introduction 



o 

. 

, Correlations of particles at low relative velocities have been used to determine sizes 

of emitting regions in heavy-ion reactions [[U [2]. In general, the smaller the emit- 
ting region the stronger are final-state effects in emission and stronger correlations 
' between the emitted particles. However, shapes of emitting regions are also of 

H \ interest. For example, a prolonged emission of particles, such as associated with 

a transitional behavior from quark-gluon plasma, can produce a relative emission 
source elongated in the direction of pair total momentum [^3J. When the final-state 
effects are due to pure identity interference, the correlation function represents a 
Fourier-transform of the relative emission source [H]. Most commonly, in the past, 
source shapes have been analyzed for charged pions only, by representing the sources 
in terms of a single anisotropic Gaussian, allowing to interpret Coulomb-corrected 
correlation functions in terms of the Fourier-transformed Gaussian [ El H] . Here, I 
discuss an alternative strategy, where the correlation function and source functions 
are expanded in surface spherical harmonics. Of the latter, the Carthesian harmon- 
ics [[S] exhibit particularly pleasing properties for the purpose. The strategy is not 
confined to correlations dominated by identity interference. 
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2. Correlations at Low Relative Velocity 

Possibility of learning on the geometry of emitting regions in reactions relies on 
the ability to factorize the amplitude for the reaction into a wavefunction $q^\ 
for the pair of detected particles, and an amplitude remnant. At at low relative 
velocity, |$q^^p may exhibit pronounced spatial features such as associated with 
identity interference, resonances or Coulomb repulsion. The wavefunction features 
can be regulated by changing the relative particle momentum q. The amplitude 
remnant, squared within a cross section and summed over the unobserved particles 
and integrated over their momenta, yields S' which represents general features of 
reaction geometry, without a significant variation with q. When examining the 
inclusive two-particle cross-section then, the structures in can be used to 

explore the geometry in S': 



da 



dpi dp2 



J dr5^(r)|<I>(-)(r)p. (1) 



The naive expectation, met at large q for a multi-particle final state, is that 
the emission of two particles is uncorrelated. It is then interesting to normalize the 
two-particle cross section with a product of single-particle cross sections and to look 
for the evidence of correlations: 

^(q) = itS^ - 1 = /dr (|<-)(r)P - l) 5p(r) . (2) 

(T dpi (7 dp2 

The second equality follows from the fact that, at large q, the relative wavefunction 
squared is equal, on the average, to 1; in combining this with 7?. ~ at large q, 
the source S, following the cross-section normalization, turns out to be normalized 
to 1, / dr S'(r) = 1. Equation ^ links the correlation function TZ to the interplay of 
deviations of the relative wavefunction from 1 with the geometry in S. Normalized 
to 1, iS* may be interpreted as the probability density for emitting the two particles 
at the relative separation r within the particle center of mass. The emission is 
integrated over time, as far as S is concerned. 



3. Correlations in Terms of Cartesian Harmonics 

Given Ti{q) and one can try to learn about S; mathematically, this represents a 
difficult problem involving an inversion of the integral kernel K{q, r) = |$q \r)\'^ — 
1, in three dimensions. The situation gets simplified by the fact that K depends 
only on g, r and the angle 9qr- Correspondingly, the kernel can be expanded in 
Legendre polynomials: 



K{(i,r)=Y,i^e+'^)Ki{q,r)P'{cose) 
I 



(3) 
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Upon expanding the correlation and source functions in spherical harmonics Y^"^, 
7^(q) = V47rX;£„7^'^™(g) Y^'"(q), one finds that the three-dimensional relation ^ 
is equivalent to a set of one-dimensional relations for the corresponding harmonic 
coefficients [|4], 

7^^™(q) = 47ry drr^ Ke{q,r) S'^'"(r) . (4) 

For weak anisotropics, only \ow-i coefficients of 7?. or S" are expected to be significant. 
The i = version of (j4]) connects the angle-averaged functions. 

The above results suggest an analysis of the correlation functions and source 
functions in terms of the harmonic coefRcients. An issue, however, is that it is 
cumbersome to analyze real functions in terms of complex coefficients that lack a 
clear interpretation for £ > 1. It becomes then natural to seek another basis for the 
directional decomposition of correlation functions and sources, one that would be 
real and have a a clear geometric meaning. Such a basis may be constructed starting 
from a unit direction vector fia = (sin0cos0, sin sine/), cos 0). A convenient choice 
of axes relative to which the angles are determined in a reaction is z along the beam 
axis, X along the transverse component of pair total momentum and y perpendicular 
to the two other axes. 

The tensor product of £ vectors fia yields a symmetric rank-^ Cartesian tensor 
that is a combination of spherical tensors of rank £' < £ and the same evenness as £: 

{n )ai...ae = f^ai f^ai ■ ■ ■ flae — ^ ^ Q'm Y . (5) 

£'<e,m 

A projection operator V may be constructed in the space of symmetric rank-^ 
tensors, out of a combination of Kronecker (5-symbols, that makes a symmetric 
tensor traceless, 

^(7'nOaaa3...a, =0. (6) 

a 

The tracelessness of (Vn^) ensures [0 1^ that the products {'Pn^)ai...at are so- 
lutions of the Laplace equation and, thus, the {'P'h^)ai...ae = Aai...ai = -^^^e^ yiy 
are combinations of spherical harmonics of rank £ only. In the above, £x, £y and £z 
are the number of times of repeated, respectively, x, y and z indices in the symmet- 
ric tensor Vh^ = A^^\ £x + £y + £z — £■ The Cartesian tensor components of A'^^'^ 
are real and may be used to replace F^™. The lowest-rank tensors are: 

-"-7 ■^a I — "ct I -^Qi Q2'' — "^1 "^2 3 ' 
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The completeness relation in terms of the Cartesian components is [[3] 

I OL\...OLi 

I a.\...OLi 

where the second equality follows from T) = = T)^ . The completeness relation 
can be utilized for expanding 7?. or S' in terms of the Cartesian tensor components 

7^(q) = 1 df7' ~ n) 7^(q') ^Y. H ■ • ■ > (9) 

£ Qi...a^ 

where the coefficients are angular moments, 

= ^/^^«(.)^S....,- (10) 

With @ and with the TZ and S Cartesian coefficients being identical combi- 
nations of the respective F^™-coefficients, the corresponding TZ and S Cartesian 
coefRcients are directly related to each other, 

n^Il...o.M) = J drr^Kiiq,r)si'l..c.,{r). (11) 

For weak anisotropies, only low-£ Cartesian coefficients for the source or correlation 

matter, 

TZ{q) = 7e(") {q) + ^ 7^(:) (g) + ^ ^"'U (?) 9a, 9.. + ■ • ■ • (12) 

The correlation coefficients T?.^^-* allow to summarize three-dimensional informa- 
tion accumulated in the measurements, in terms of one-dimensional plots. With Q 
and p2p . the values of correlation (and similarly source) function in the +z direction 
are given by 

7^(g) = 7^(°)(g)^-7^«(g) + 7^i2)(g) + ... 

and in the xz plane at 45° degrees by 

niq) = 7^(")(g) + ^7^(l)(g) + ni'^q) + \ {n^^KD + T^fK<D) + ■■■■ 

The coefficients TZ'^^^ can be found through angular integration ([TU| or by fitting 
the expansion series ([12]) to the dependence of data on spherical angle. The latter 
procedure can be particularly useful for an experiment blind in certain directions. 
For relating the correlation coefficients to source coefficients, the two-body problem 
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needs to be solved for the wavefunction $q -^(r). Not only identity interference, but 

also Coulomb and strong interactions give rise to structures in |$q~''(r)p that yield 
finite kernels Ki for multipolarities £ > [[5]. Contrary to naive expectations, even 
strong interactions acting in s-wave only produce finite kernels a.i I > 0. With 
regard to the Coulomb interactions, there is no good reason trying to correct for 
them [[2], as they can provide access to source deformations just as the identity 
interference can. Features of the source can be accessed in pT|) by assuming a 
parameterized shape of the source and fitting source parameters to the correlation 
coefficients, or by imaging source coefficients in the same manner as for £ = [[7]. 

4. Sample Results from Data Analysis 

The number of angular expansion coefhcients that need to be considered for the 
correlation and source functions can be quite small. Thus, for pairs of identical 
particles, such as identical pions, the correlation functions are symmetric under 
inversion of relative momentum, q — > —q, meaning that only even-^ components 
are finite. Around midrapidity of a symmetric system, such as Pb + Pb or Au + 
Au, the correlation function should be symmetric with respect to the interchange of 
forward and backward directions, qz —> —Qz, meaning that only even-z components 
of the functions can be finite. When the correlation functions are averaged over 
reaction-plane orientation, they are symmetric under interchange of the orientation 
of sideway axis, qy —^ —Qy, meaning that only even-y components of the functions 
can be finite. In the end, as a consequence of the above, the correlation function in 
the considered situation must be invariant under the interchange of the orientation 
of outward axis, qx — > —qx, meaning that only even- a; components of the functions 
can be finite. For £ < 4, the finite Cartesian moments are then £ = 0, x'^ , y'^ and 
for £ = 2 and x'^, y*, z^, y^, x^ and y^ for £ = 4. Due to the tracelessness 
of the tensors from Cartesian components, however, only 2 of the above t = 2 
components are independent and only 3 of the £ — A. Thus, depending on data 
statistics, the dependence of anisotropic correlation function on spherical angle at 
any q can be described with either 3, for ^ < 4, or 6, for ^ < 6, functions of relative 
momentum q. 

Figure [1] shows 1 = and 1 = 2 Cartesian harmonic results for the midrapidity 
TT^-TT^ correlation functions measured by the NA49 Collaboration in central Pb -|- 
Pb collisions at ^ = 17.3 AGeV. The £ = panels show the correlation function 
averaged over directions of relative momentum and the i = 2 panels illustrate the 
quadrupole deformation of the function. The coefficients are given by negative of 
the sum of and ?/2 coefficients, 7^^7 = -{T^x2 +7^^7)■ No Coulomb corrections 
have been applied to the data. On the other hand, effects of Coulomb interactions 
are accounted for directly in relating source functions to the data. In the left panels 
of Fig. [TJ it is seen that a single Gaussian source fails to describe the data at low 
relative momenta. In the right panels, on the other hand, it is seen that a source 
combining two Gaussian sources provides a much better description of the data. 



6 



P. Danielewicz 



NA49 PbPb vs=±f7.3 AGeV 



0.2 



X 0.1 



^™ 0.05 

or 



of 



0.05 



• Data Kj „ 

Ellipsoid Fit 0<cen<7% 
ftSg 0<Py<70MeV 

^'■»iB ■0.35<yL-yo<0.35 
.'•f».««»»««»«»»»»«»t 



b) 



NA49 PRELIMINARY 



jj?.*.*?"**fl*Si 



c) 



A = 0.200 ±0.003 
r,^= 5.34± 0.04 fm 
r = 4.92±0.04fm 
= 7.22± 0.07 fm 



:4.9 



NA49 PbPb \s=±17.3 AGeV 



0.2 



a: 0.1 



- 



0.05 



- 



0.05 

DC 



- 



0.01 0.02 0.03 0.04 0.05 
q (GeV) 



0.06 



• Data KK 
LI 2-SrcFit 



0<oer\<7% 
l«a o<PT<70MeV 
■0.35<yL-yo<0.35 





b) NA49PREUr,^f^|y^QQ^ 

r' = 3.8±0.l'fm 
"f r'= 4.2±0.1 fm 



, r;= 3.7±0.2fm 



■e 



= 0.32 ± 0.01 
r, = 8.1±0.2fm 
r = 6.5+ 0.2 fm 



yr 



I 



12.6± 0.4 fm 
7ndf = 1 .2 



0.01 



0.02 0.03 0.04 
q (GeV) 



0.05 0.06 



Fig. 1.^ = and two independent £ — 2 angular moments of tt^-tt^ correlation 
functions in central Pb + Pb collisions at ^/s = 17.3 A GeV. Filled circles represent 
measurements of the NA49 Collaboration [ [8] . Squares in the left panels represent 
a fit to the data with a source in the form of an anisotropic Gaussian. Squares in 
the right panels represent a fit to the data with a source combining two anisotropic 
Gaussians. 



Figure [2] shows next source values along x, y and z axes, from summing £ = 
and £ = 2 source coefficients deduced from the coefficients of correlation function 
measured in central Pb + Pb collisions at y/s — 8.9 AGeV and ^/s — 17.3 AGeV 
[[8]. Aside from the single and two Gaussian results, also source values obtained 
through imaging [i7j are shown. The imaged sources usually yield a very good de- 
scription of measured correlations. It is seen in the figured that the fitted 2-Gaussian 
and imaged sources agree fairly well with each other. At both energies, those sources 
exhibit tails extending over 30 fm of relative distance, in the longitudinal z-direction. 
The failure of the single-Gaussian source to reproduce the correlation moments at 
low g, observed in Fig.[Tl appears to be related to the failure in producing such a tail 
at large r. A non-Gaussian enhancement in the outward a;-direction further devel- 
ops in the imaged and 2-Gaussian sources between the energies of = 8.9 AGeV 
and y/s = 17.3 AGeV. The sources in Fig.[2]are most compact in the y direction. A 
tail in the source in the z direction is expected in connection with a strong collective 
expansion of the system along the beam [ [6] . Greater extension in the x- than the 
y- direction and the tail in the x-direction can be associated with emission extended 
over time. Low transverse momenta have been purposely selected for the analysis 
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to eliminate any significant effects of Lorentz 7-factors. 



NA49 Pb+Pb •is=ta.9AGeV NA49 Pb+Pb \ s=± 17.3AGeV 




r (fm) r (fm) 



Fig. 2. Relative midrapidity tt^-tt^ source in central Pb + Pb collisions along 
X, y and z coordinate axes, from summing the £ < 2 angular components of the 
source. The left and right panels represent, respectively, source in collisions at 
-^s — 8.9 AGeV and ^/s = 17.3 AGeV. Squares represent source values obtained 
by imaging the angular components of correlation functions measured by the NA49 
Collaboration [^SJ. Triangles represent results of the fit to the measured correlation 
moments [[9] with a source in the form of an anisotropic Gaussian. Circles represent 
results of the fit to the measured correlation moments with a source combining two 
anisotropic Gaussians. 

Extended non-Gaussian tails in the longitudinal z-direction are further observed 
in central collisions at RHIC energies. Figure [3] shows the source values extracted 
from the analysis of £ < 4 angular moments of midrapidity correlations measured 
by the PHENIX Collaboration in central Au + Au coUisions at a/s = 200 A GeV. 
Again, a fit with a single Gaussian to the correlation functions fails to produce the 
tail along z-axis emerging from the imaging. The imaged source also exhibits some 
tail in the outward direction which is not there in the single Gaussian source. 

Figure [3] compares further the sources extracted from PHENIX data to those 
predicted by model calculations. The Therminator model [HO], represented in the 
left panels, has been successful in the past to reproduce Gaussian-source parameters 
extracted from correlation sources corrected for Coulomb effects. The model sup- 
plements short-range features of a source with long-range features due to resonance 
decays. In the left panels of Fig. [3l it is observed that the model can fully explain 
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only the experimental source values in the y-direction and partly in the z-direction. 
However, the model fails completely as far as the long-range features of the imaged 
source in the x direction. The latter is quite independent of whether the resonance 
decays take place or are suppressed. 

The multi-phase transport model (AMPT) [ [11] follows the dynamics of both 
partonic and hadronic phases and includes resonance decays. As is apparent in the 
right panels of Fig. [3l the model describes fairly well the imaged source features 
along the x and z directions but fails rather dramatically with respect to the y 
direction. 



PHENIX Run4 Au+Au ^s=200AGeV PHENIX Run4 Au+Au ^s=200AGeV 




r (fm) r (fm) 



Fig. 3. Relative midrapidity 7r+-7r+ source in central Au -I- Au collisions at 
y/s = 200 A GeV, along x, y and z coordinate axes, from summing the £ < 4 
angular components of the source. Squares represent source values obtained by 
imaging the angular components of correlation functions measured by the PHENIX 
Collaboration [[9]. Triangles represent results of the fit to the measured correlation 
moments [[5] with a source in the form of an anisotropic Gaussian. In the left pan- 
els, circles represent a fit using the Therminator model [UHl [H] with, respectively, 
resonance decays operating, for full circles, and switched off, for open circles. In 
the right panels, circles represent results of the AMPT model [lllj. 



5. Conclusions 

As demonstrated with examples of the analyses of NA49 and PHENIX correlation 
measurements, the expansion in Cartesian harmonics represents a practical strategy 
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for summarizing correlation functions and for accessing sources shapes. The strategy 
is capable of producing a wealth of detail on emission sources and challenging this 
way theoretical interpretations of the data. 
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